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ABSTRACT 

The presence of anthropogenic RFI is expected to adversely 
impact soil moisture measurement by NASA’s Soil Moisture 
Active Passive mission. The digital signal processing ap- 
proach and preliminary design for detecting and mitigating 
this RFI is presented in this paper. This approach is largely 
based upon the work of Johnson [1] and Ruf [2]. 

Index Terms — SMAP, Soil Moisture Active Passive, 
Digital Radiometer, Filterbank, Kurtosis 

1. INTRODUCTION 

Radio-frequency interference (RFI) has detrimentally af- 
fected soil moisture remote sensing data since 1978 when it 
was first observed at 6.6 GHz with SMMR on SeaSat [3]. 
Over the last several years it has started to seriously impact 
remote sensing science products and is expected to con- 
tinue. This is primarily because the frequency requirements 
of remote sensing have begun to overlap active radio service 
allocations. Until recently, soil moisture was being sensed in 
the C-band, for which there is no passive service allocation. 
The launch of the European Space Agency’s Soil Moisture 
Ocean Salinity (SMOS) mission [4] begins a new era of L- 
band remote sensing from space, which has been long desired 
for remote sensing of soil moisture through vegetation. The 
L-band record from space will be continued with NASA’s 
Aquarius [5] and SMAP [6] missions. 

Although there is a primary exclusive allocation to pas- 
sive use in the L-band, it is still highly susceptible to RFI 
because of the proliferation of the active services in border- 
ing frequency allocations. Thus, RFI detection and mitigation 
strategies are needed so SMAP can meet its requirements. 

The chosen approach for SMAP is to use a digital process- 
ing backend in lieu of microwave power detectors. Johnson 
produced the first RFI-mitigating digital radiometer [1] and 
Ruf, et al., introduced the kurtosis statistic as an RFI detec- 
tor and classifier [2]. Using their previous research plus work 


done for the Hydros mission [7] as its basis, the SMAP dig- 
ital backend will provide both time and frequency diversity 
data at each footprint to enable the ground processing algo- 
rithms to detect and remove RFI. To do so, it will use digital 
signal processing (DSP) techniques to measure 1200 samples 
in time and frequency for each SMAP footprint. This paper 
summarizes the DSP algorithms and electronics implemen- 
tation using field-programmable gate arrays (FPGAs) in the 
SMAP preliminary design. 

2. SMAP HARDWARE IMPLEMENTATION 

2.1. Radiometer Front End (RFE) Analog Subsystem 

The SMAP mission uses a single feed polarimetric radiometer 
that operates in the protected L-band radio frequency range 
from 1400 to 1427 MHz. The received radiometer signal, 
x(t ), is branched into horizontal and vertical polarizations, 
denoted by E H (t) and E v (t ), respectively. In each branch, 
the polarized signals are bandlimited, amplified, and down- 
converted to a 120 MHz intermediate frequency (IF). The 
Radiometer Digital Electronics (RDE) subsystem digitizes 
Eh (t) and E v ( t ), forming digital signals E H (n) and E v (n) 
respectively. At the core of the RDE is a FPGA-based DSP 
system that computes time-frequency diversity data used for 
subsequent ground processing, removal of RFI, and geophys- 
ical parameter retrieval. 

2.2. Digital Signal Processing Approach 

SMAP will be the first spacebome mission to use on-board 
DSP for RFI mitigation using state-of-the-art algorithms de- 
veloped by Ruf [1] and Johnson [2]. In prior spacebome 
radiometers, the IF signal power is measured by an analog 
square-law detector followed by a lowpass filter. The result- 
ing power measurement is digitized and transmitted to the 
ground. The SMAP radiometer digitizes the IF signal di- 
rectly and uses DSP to compute IF signal power, as well as 
additional data products useful for RFI time-frequency local- 
ization and statistical detection using kurtosis. 



The SMAP DSP computes the first four raw sample mo- 
ments of the fullband signals, E H (n) and E H (n) . In addition, 
the fullband signal is split into 16 adjacent and overlapping 
subband signals ejy (n) and (n). The first four raw sam- 
ple moments are computed by the DSP for each subband sig- 
nal as well. Finally, the cross-correlation of the fullband and 
subband signals are computed. The integration times used for 
each fullband and subband data product are programmable, 
but are initialized to a minimum of 280 and 1120 /is, respec- 
tively. All of these statistics are transmitted to the ground and 
are used as inputs to the RFI detection algorithm. The system 
functional block diagram is shown in Figure 1. 

2.3. Radiometer Digital Electronics Subsystem 

The SMAP RDE consists of three processor boards. The 
first two identical boards are called Analog Processing Units 
(APU)s. These boards, named APU-H and APU-V, digitize 
the input signals with 14 bits at a sample rate of 96 MHz. The 
third processor board is called the Digital Processing Unit 
(DPU). The DPU collects and prepares data products from 
each APU card for ground transmission. The DPU also com- 
putes cross correlation between horizontal and vertical polar- 
ized fullband and subband signals. 

2.4. FPGA-based DSP 

On-board radiometer DSP is performed using a total of seven 
antifuse-based Actel RTAX2000 FPGAs - three on each APU 
board and one on the DPU board. If SRAM-based FPGAs 
with significantly more logic gates are used, the total number 
of FPGAs can be reduced to about 1 or 2 devices. However, 
factors such as flight qualification, radiation hardness and 
flight heritage deemed these kinds of FPGAs too risky for 
spacebome application. 

Actel RTAX4000 FPGAs, which have about twice the 
number of logic gates as the RTAX2000, offered a radiation 
hardened solution, but the risk of using these parts was also 
determined to be unacceptably high. This is because they had 
no prior flight heritage; had no easy and low-cost prototyping 
solution; and were too early in the flight qualification stage to 
be considered viable. 

With the constraint of using Actel RTAX2000 devices 
to do the signal processing, several algorithm optimizations 
were employed that allowed for a logic-area efficient imple- 
mentation of the RFI mitigation algorithm. 

A number of DSP techniques [8, 9, 10] commonly found 
in digital communication receivers and radar are applied to 
the radiometer DSP algorithm to minimize FPGA logic area. 
Algorithms amenable to FPGA multiplier sharing and signal 
sample interleaving were deliberately sought out and are used 


extensively throughout the entire processing chain to mini- 
mize FPGA logic. 

3. PROCESSING SUBSYSTEMS 

The DSP is divided into five modules: Signal Pre-Processor 
(SPP), Digital Down Converter (DDC), Polyphase Filter 
Bank (PFB), Statistics Calculation Unit (SCU) and Cross 
Correlator (XCORR). Each module discussed in the fol- 
lowing sub- sections. In the follwing discussion, denote the 
fullband signal (from APU-H or APU-V) by E(n) and the 
subband signal by (n). Subscripts H and V are used only 
when making a distinction between the polarized signals. 

3.1. Precision and numeric representation 

All FPGA computations are preformed using signed, 2’s- 
complement fixed-point arithmetic. Full-precision fixed point 
is maintained in each processing module, but the word length 
is limited to 20 bits between modules. A finite-precision 
rounding system similar to block floating point is employed 
throughout the DSP, where a variable-start, contiguous range 
of 20 bits is selected from the full precision bus with the aid of 
in-orbit programmable starting pointers. Due to limited trans- 
mit bandwidth to the ground, each data product generated by 
the RDE is restricted to 12 bits. 

3.2. Signal Preprocessor 

It is desirable to have robust diagnostic capabilities for the 
RDE once in orbit. Two diagnostic capabilities in the DSP 
are managed by the SPP. First, the SPP allows the user to 
bypass the radiometer signal and inject test signals instead. 
The SPP also permits the user to add a programmable DC- 
offset value to the input signal to compensate for DC bias. 
A third diagnostic feature of the RDE bypasses processing 
all together, allowing short time records of E(n) to be saved 
directly and transmitted to the ground for analysis. 

3.3. Digital Downconverter 

The RDE receives a 24 MHz bandlimited analog signal cen- 
tered at 120 MHz IF. The signal is subsampled at 96 MHz, 
allowing the alias in first Nyquist zone to be digitized. The 
DDC downconverts E(n ), reducing its sample rate to 24 
MHz. The mixer stage, image rejection filter, and down- 
sampler are all combined into an efficient periodically time- 
varying (PTV) [11] multirate FIR filter using an approach 
similar to Ellingson [12]. This architecture minimizes FPGA 
logic utilization by reducing the number of required multipli- 
ers from 130 to 16. The output of the DDC is referred to as 
the fullband signal denoted by E{n). 

Although the RFE aggressively bandlimits the L-band sig- 
nal to 24 MHz using two 8-pole filters, it is important to have 
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Fig. 1. SMAP-RDE signal processing flow. Output data consists of fullband and subband moments for each polarization and 
each I/Q component, as well as fullband and subband cross correlation. 


high quality image rejection in the DDC. The image rejec- 
tion filter for the DDC is a 64-tap 12-bit coefficient FIR low- 
pass filter based on a Blackman window with -3dB point at 
12MHz. This filter suppresses the mixer image and out-of- 
band interference down to -60dB. 


3.4. Critically-Sampled Polyphase Filter Bank 

The PFB splits E(n) into 16 uniformly spaced, equal- 
bandwidth subband signals e^(n), k = 0, ...,15. The 
filter bank requirements are listed in Table 1, leading to the 
response in Figure 2. The PFB design requirements are met 


Table 1. SMAP Filter Bank Specifications 


Channels 

16 

Total bandwidth 

24MHz 

B andwidth/channel 

1.5 MHz 

k to k + 2 th stopband attenuation 

32 dB 


with a prototype 64^-order unity gain FTR filter based on 
a Kaiser window with /3 = 4 and -2.8dB overlap point of 
750 kHz. The coefficients are quantized to 12-bits to achieve 
about 60 dB stopband attenuation. 

The filter bank is comprised of two units. The first is a 
PTV FIR filter that cycles through the coefficients of each 
polyphase sub-filter. The second is a standard 16-point, 
pipelined FFT. Compared with a uniform filterbank and a 
parallel-branch implementation of a polyphase filter bank [9], 
this filter bank architecture minimizes FPGA logic. 

Much like the DDC, the desired side effect of this archi- 
tecture is the interleaving of subband signals. Every block of 
16 serialized samples represents a time sample from each of 
the 16 subbands. The output of the filter bank is denoted as 
e t fc ](n) = i^(n) + jq^ k \n), k = 0, . . . , 15. 


Kaiser (p=4), Window Length = 64, 12-bit coefficients 
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Fig. 2. Frequency response of polyphase filter bank 

3.5. Statistics Calculation Unit 

The role of the SCU is to compute the first four raw sample 
moments of E(n) and e^(n). The r th moment, m^\ r = 
1, . . . , 4 for the in-phase and quadrature component signals 
for each subband are computed in Equation 1, 

= av,fc ^2 (* [fe] (n)) (1) 

71=0 

where a accounts for scaling the full-precision moment com- 
putations down to fit within a 12-bit bus subset. Similarly, the 
moments M r of the E(n) are computed using an equation of 
the same form as above by dropping the superscript [k] and 
using capital notation. 

3.6. Cross Correlator 

The last DSP module in the SMAP-RDE signal processing 
chain is the cross-correlator. Cross-correlation necessarily 














takes place on the DPU board. The subband cross correla- 
tion between e ^ (n) and (n) is given by, 

r HV = ( e H ] ( n ) (fiy'w) ) 

= 7 (*H ( n )*F ( n ) + 9 h ( n )«v (»)) 

71=0 

+ J7 X] (*v ( n ) “ ( n )«y (»)) . ( 2 ) 

71=0 

where k = 0, . . . , 15, and 7 represents a scaling factor intro- 
duced by selecting a 12 -bit subset of the full precision fixed 
point result. The complex cross-correlation of the fullband 
signals E H (n) and E v (n) are computed with and expression 
similar to Equation 2 by dropping the superscript [k] and us- 
ing upper-case notation for the signals. 

4 GROUND ALGORITHMS 

Using the detailed information provided by the SMAP in- 
strument, the ground processing algorithms will detect RJFT 
in the raw data, remove it, and output calibrated nearly RFI- 
free brightness temperature estimates. The diversity of the 
SMAP data product allows several different detectors to be 
utilized. A time domain pulse detector can be used on the 
full band product for removing radar interference. A cross- 
frequency detector can be used to remove persistent narrow 
band sources. It is envisioned the kurtosis products will aid in 
classifying RFI, if not in detecting less- well-behaved sources. 
A maximum likelihood decision process will be used to com- 
bine the detector outputs and produce a single estimate of an- 
tenna temperature with RFI removed. 

5. SUMMARY 

The development of the SMAP RFI-mitigating radiometer is 
based on the ground work of Johnson and Ruf. Application of 
fundamental DSP tools, novel algorithm simplifications and 
FPGA architecture tricks make this digital radiometer possi- 
ble for space implementation. The preliminary design proto- 
type is being tested currently using an Altera Stratix m DSP 
kit in preparation for a preliminary design review in Septem- 
ber 2010. The spaceflight engineering test unit, currently be- 
ing designed, is scheduled for delivery in 2011 . 
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